Anderson localization in optical lattices with speckle disorder 
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We study the localization properties of non-interacting waves propagating in a speckle-like po- 
tential superposed on a one-dimensional lattice. Using a decimation/renormalization procedure, 
we estimate the localization length for a tight-binding Hamiltonian where site-energies are square- 
sinc-correlated random variables. By decreasing the width of the correlation function, the disorder 
patterns approaches a 5-correlated disorder, and the localization length becomes almost energy- 
independent in the strong disorder limit. We show that this regime can be reached for a size of the 
speckle grains of the order of (lower than) four lattice steps. 
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I. INTRODUCTION 



The Anderson model was proposed first to explain the 
absence of electronic diffusion in certain random lattices 
[l[ and later to explain the absence of diffusion of light 
in certain amorphous materials [2j. In the first case, 
the electron energy is lower than the maxima of the lat- 
tice potential and the particles diffuse by tunnelling. In 
the second case the particles (photons in Q) energy is 
higher than the potential and the particles are "free". 
In both cases, in the absence of disorder, the eigen- 
states are delocalized states (Bloch waves in the first 
regime, plane waves in the second). The presence of a 
(^-correlated random potential freezes the wave propaga- 
tion at a length of Ci oc , the localization length, always 
in one (ID) and two dimension (2D), and depending on 
the disorder strength and the energy wave also in three 
dimension (3D) [3J. In the last few years the ultracold 
atom community has devoted a large effort to the exper- 
imental realization of Anderson localization. Anderson 
localization in the tight-binding regime was observed in 
momentum space with kicked- rotor set-ups, in ID 
and 3D [5J, and in real space by using a quasi-periodic 
potential (thus not strictly-speaking a random potential) . 
The key-ingredient for the experimental study of Ander- 
son localization of ultracold atoms in the "free" -particle 
regime has been the speckle potential [fj. By using this 
handable optical potential, Anderson localization was ob- 
served in ID [7j and in 3D |8j,|9|, and anomalous diffusion 
was observed in 2D [10( ■ The auto-correlation function of 
the speckle potential is a square-sinc, thus it decays alge- 
braically as a long-range correlated disorder, but is also 
characterized by a finite size correlation length w, corre- 
sponding to the width of central bump of the square-sinc 
function. 

In ID, the presence of disorder patterns with auto- 
correlation functions decaying algebraically can mimic 
the presence of a mobility edge [Ill4l4{ and can even en- 
hance localization, as shown in a microwave experiment 
1 1( | . Both phenomena are due to the fact that the disor- 



der spectrum is non-zero in a finite momentum interval. 

Very recently Semmler and coworkers [l5j have stud- 
ied the phase diagram of correlated fermions in 2D and 
3D optical lattices and in the presence of a speckle po- 
tential. From the analysis of the local Density Of States 
(DOS) they identify an Anderson-Mott and a Mott lo- 
calized phase as functions of the interaction strength 
and the strength of the speckle potential. In this arti- 
cle we analyze the possibility of observing Anderson lo- 
calization of a non-interacting wave, for example a non- 
interacting Bose-Einstein condensate [l6|, in a speckle 
potential superposed to a ID lattice potential. By us- 
ing a decimation/renormalization scheme [l7j we analyze 
how the DOS of a lattice is modified by the presence of 
the speckle, and we estimate the localization length Ci oc 
as function of the disorder strength and of the width w 
of the auto-correlation function. The speckle potential 
is introduced as an on-site disorder which has statistical 
properties which are the same as a genuine speckle poten- 
tial. This is illustrated in Sec. [TTJ In Sec. IIIII we remind 
the reader of the decimation/renormalization procedure 
exploited to compute the DOS and Ci oc . Our results 
show how the efficacy of the speckle potential to localize 
increases by increasing the disorder strength and by de- 
creasing the correlation length w. These results can be 
a guide to choosing the experimental parameters to ob- 
serve Anderson localization in the tight binding regime 
with speckle disordered patterns. 



II. THE MODEL 

To study the effect of a speckle potential in the pres- 
ence of a ID lattice on matter- wave transport we use the 
ID Tight-Binding (TB) Hamiltonian, 

n s n s — 1 

H = J2E i \i)(i\+ ^ t (\i)(i+l\ + \i+l)(i\) (1) 

i=l i=l 

where n s is the number of sites, Ei the energy at the site 
i. The hopping term t is chosen site- independent. The 
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effect of the speckle potential is introduced in the on-site 
energy distribution by setting 



Ci = (SEiSE i+e ) = s 2 



sm(2n£/w) 
2-Kijw 



where SEi = Ei — (Ei) is the fluctuation of Ei with re- 
spect to the mean value (Ei), s = yj ((SEi) 2 ) is the disor- 
der strength and w the width, in the units of the lattice 
step d, of the correlation function in Eq. @. The dis- 
order spectrum is not uniform as in the Anderson model 
[l|, but is described by the triangular function Q 

S k oc s 2 (re- \k\)6(n - \k\) (3) 

where re = 4tt/w, and 0(x) is the Heaviside function. 

A. Generation of the disordered potential 



We use the Fourier Filtering Method (FFM) [lj-|2j| 
to generate the disorder pattern described by the cor- 
relation function @. First we generate a sequence of 
N 5-correlated random numbers {uj}, with j = 1, . . . N 
from a uniform distribution centered in zero and of width 
1. The second step is the generation of the desired {Ej} 
distribution by "filtering" , in Fourier space, the uniform 
distribution {uj}. The filter being the spectral function 
Sk , the Ej 's are evaluated directly from the expression 



N k N 



E; 



FEE VS~ k e^h 



(4) 



jk=0 m=l 



where Nk = 8N/w et k = — k+(tt /N)jk- By construction 
(Ej) = 0, namely SEj = Ej, and (EjEj + g) verifies Eq. 
([2]) in the limit N — > oo. A different choice of (Ej) would 
just shift the zero of the energy. 



III. NUMERICAL RESULTS: THE DOS AND 
THE LOCALIZATION LENGTH 

In the continuous limit, the single-particle DOS for an 
optical speckle has been studied in [2l| . In the presence of 
a lattice (and in the absence of the speckle potential), the 
low-energy single-particle DOS has a typical saddle shape 
with two horns that correspond respectively to the center 
and the edge of the first Brillouin zone. To evaluate how 
the speckle potential modifies the DOS of the lattice, we 
compute the DOS, M(E), regarding the Hamiltonian ^ 
by using the Kirkman-Pendry relation [22j 



Af(E) 



lim — Im 



1 f d ln[Gi, tt . (E + ia)} 



dE 



(5) 



Here G(E) = (E — H) 1 is the Green's function re- 
lated to the Hamiltonian H at energy E, and Gij(E) = 



(i\G(E)\j). With the aim of computing the matrix ele- 
ment Gi^n a (E), we reduce the dimensionality of the sys- 
tem by evaluating the effective Hamiltonian 



(2) H = E 1 \l)(l \+E ns | n s )(n s \+t(\ l)(n s + l \ + \ n s + l)(l |), 

(6) 



where E\, E ns and t are functions of the energy E and 
of the Hamiltonian elements of the decimated states (2, 
3, . . . n s - 1) (HI S3. The Green's function of the effec- 
tive Hamiltonian ©, G(E) = (E — H)" 1 , coincides with 
G(E) in the subspace {l,n s } by construction. 

The numerical results for the DOS as a function of the 
energy in units of \t\ are shown in the first column of 
Fig. [TJ One can observe that for large values of w, the 
speckle disorder mainly affects the edge states of the DOS 
of the underlying perfect chain, while for smaller values 
of w the disorder mainly influences the central part of 
the spectrum. 

The presence of the disorder modifies not only the 
DOS, but also the nature of the states, from extended 
to localized. In the continuous limit, the presence of 
the correlations described by Eq. ^ does not destroy 
localization but deeply modifies the behaviour of the lo- 
calization length as a function of the energy [l2T4l4l |. To 
study the behaviour of the localization length Ci oc (E) 
in the tight-binding regime, we compute the Lyapunov 
coefficient j(E), through the asymptotic relation 



1 



>y(E) = lCioc(E)}- 1 = lim —In 

n s ->oo n s U 



G r , 



,(E) 



Gl,n 3 (E) 



(7) 



The results shown in the second column of Fig. Q] have 
been computed for the case n s — 200, but we have 
checked that the values obtained do not change signif- 
icantly by increasing the value of n s up to 1000. Analo- 
gously to the continuous case, we observe that all states 
are localized. In the limit of weak disorder, the localiza- 
tion length at the center of the spectrum, Ci oc (E = 0), 
is quite large, of the order of 50 lattice sites. By in- 
creasing the strength of the disorder, £i oc (E = 0) de- 
creases significantly only for small values of the correla- 
tion length (w — 7r and w — 2tt/3), and Ci oc (E) becomes 
almost energy-independent in the whole band. Longer- 
range correlations (w — 2tt and w — 47r) act instead more 
efficiently on the edge states. To better understand these 
reasults we can refer to the continuous case, where 



£ioc(k)~ 



A s oc (fe)" 



w 



2k- 



(8) 



in the Born approximation (see for instance (Tlj)- From 
Eq. ([8]) we can expect to observe (i) a decrease of the 
localization length for large values of w in the limit k — > 0, 
since Cf oc (k — ► 0) ~ k 2 /(s 2 w), and (ii) an increase in 
the localization length for k > re/2, namely where the 
Born approximation is no longer valid [l2l [l3j]. Since 
in the TB case there is a k — > n/d — k symmetry in 
the DOS due to the presence of the underlying lattice 
and correlations act symmetrically with respect to the 
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FIG. 1: (Color online) DOS and localization length Ci oc (in 
units of the lattice step d) as functions of the particle energy 
(in units of the energy hopping) for a lattice of 200 sites where 
site energies are Cf-distributed random variables [see Eq. ((2])]. 
Each curve correponds to an overage over 100 configurations. 
The different lines correspond to: w — 4n (continuous red 
line), 2-7T (dashed green line), ir (point-dashed blue line) and 
2n/3 (pointed magenta line). From top to bottom: s/\t\ = 
1,2, 5 and 10. 



center of the spectrum, the observation (ii) leads to the 
conclusion that, if one wants to experimentally observe 
Anderson localization in the whole low-energy band, k/2 
should be greater than 7r/2d, thus w should be lower than 
4 lattice steps. This finding, deduced from propagation 
in the continuous space, is in good agreement with the 
numerical results for the TB case. 



IV. CONCLUSIONS 



In this article we have studied the effetcs of a speckle 
potential on the spectrum of a quantum particle (or a 
non-interacting wave) in a lattice potential. At fixed, 
large disorder strength (s = 10|i|), the localization ef- 
ficacy of the speckle potential depends strongly on the 
width of the auto-correlation function w. Large val- 
ues of w enhance localization at very low energies and 
at the edge of the Brillouin zone. Shorter-range cor- 
relations (w < 4 lattice sites) act more efficiently on 
the center of the spectrum. More generally, our results 
show that a speckle superposed to an optical lattice is a 
suitable potential to study Anderson localization in the 
tight-binding regime: analogously to the continuous case, 
speckle correlations deeply modify the behaviour of the 
localization length as a function of the energy, but do not 
induce an insulator-metal transition. 
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